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Abstract. Subbarrier fusion of the 7 Li+ 12 C reaction is studied using an antisymmetrized 
molecular dynamics model (AMD) with an after burner, GEMINI. In AMD, 7 Li shows ana + t 
structure at its ground state and it is significantly deformed. Simulations are made near the 
Coulomb barrier energies, i.e., E cm = 3 — 8MeV. The total fusion cross section of the AMD + 
GEMINI calculations as a function of incident energy is compared to the experimental results 
and both are in good agreement at E cm > 3MeV. The cross section for the different residue 
channels of the AMD + GEMINI at E cm = 5M eV is also compared to the experimental results. 



1. Introduction 

Availability of radioactive beam facilities has stimulated theoretical and experimental interest 
in the structure of nuclei far from the stability line. Nuclear fusion reactions near the Coulomb 
barrier are strongly affected by the structure of the interacting nuclei, especially with weakly 
bound neutrons [TJ. Some theoretical calculations predict that the fusion cross section is 
enhanced over well-bound nuclei because of the larger spatial extent of halo nucleons [2j. On 
the other hand halo nuclei can easily break up in the field of the other nucleus, due to their 
low binding energies, before the two nuclei come close enough to fuse and carry away available 
energy. Early calculations of Hussein et al. [3j indicate that the actual fusion cross section 
decreases significantly at all energies. However recent couple channel calculations of Hagino 
et al. [4j have concluded that the fusion cross section increases below the Coulomb barrier 
because of the neutron halo whereas it decreases above the barrier because of weak coupling of 
the halo nucleons. Experimentally this is still a hot debate because of experimental difficulties. 
Another interest we propose here is the influence of the cluster structure in the fusion mechanism. 
The cluster structures of light nuclei have been studied theoretically from stable to those far 
away from the stability line. Recent calculations, using an antisymmetrized molecular dynamics 
model (AMD), indicate that light nuclei exhibit variety of distinct cluster structures El [S] . 



The cluster structures are predicted even for nuclei with Z ~ N of Li and Be [6] (where Z 
and N are the charge and neutron number in a nucleus, respectively). When nuclei with a well- 
developed cluster structure are involved in fusion reactions near the barrier, it will be reflected on 
the fusion cross section, especially in the variety of the exit channel distribution or in the fusion 
residue distribution. The cluster structure effect may be enhanced in the fusion reaction between 
light systems. In these systems the Coulomb interaction becomes small and the proximity effect 
between the two nuclei will be enhanced and therefore the structure of the projectile and/or 
target may reflect the fusion cross section directly, especially as an enhancement of particular 
incomplete fusion channels. In this paper we present the calculated results in the study of the 
fusion reactions of the 7 Li + 12 C system near the Coulomb barrier using AMD simulations. 

2. AMD simulations 

The initial nuclei of Li and 12 C were produced by the AMD code of Ono et al. in refs. jS], using 
the Gogny interaction. The binding energy and root mean square radius of these initial nuclei 
are compared to the experimental values in Table 1 . All calculated values are in good agreement 
to those of the experimental values, except for the root-mean square radius of 7 Li in which 
the calculated value is about 20% larger than that of the experiment. The more sophisticated 
calculation in ref. [6] of the experimental root mean square radius of 7 Li is also well reproduced 
with a distinct a + t structure. 



Table 1. initial nuclei. 
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Figure 1. Li initial nucleus with 
a + t structure. 
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Figure 2. C initial nucleus with 
3a structure. 



The density plot of these nuclei are also shown in Fig{T] and Figj2j Symbols indicate the 
location of all nucleons. One can see in both figures that nucleons are well clusterized in space. 



In the Li case, two clusters are observed, the larger one corresponds to an a and the other to 
a triton, and the nucleus is very deformed. In 12 C, 3a clusters are observed, but the nucleus is 
compact and much more spherical. 

Using these initial nuclei, 7 Li + 12 C reactions were simulated at center of mass energies 
between 3 to 8 MeV. Calculations were performed in the impact parameter range, b, from to 
7fm. In b > 7fm, no fusion reactions are observed. In Figj3] the time evolution of the density 
distributions is shown as typical examples of the complete and incomplete fusion reactions. On 
the left panel, a complete fusion reaction is observed. In the middle, only the a particle is 
transferred into the 12 C nucleus and triton is escaped as a spectator. On the right panel, only 
triton is absorbed and the a particle becomes a spectator. The latter two cases are mainly 
observed at larger impact parameters. In each incident energy, a few thousand to ten thousand 
events are generated, depending on the fusion cross section, proportional to the impact parameter 
in the given range. 
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Figure 3. Time evolution of the 2D density plots for typical fusion reactions. Impact 
parameter, incident energy and reaction product at the bottom of the simulations are indicated 
on the top of each figure. The density plot is made by projecting that of all nucleons on the 
X-Z plane. The contour lines are plotted on a linear scale. 

The AMD calculations were performed up to times ranging from 3000 fm/c at lower energies to 
1000 fm/c at higher energy side and clusterized at the end of the calculation, using a coalescence 
technique in phase space. The coalescence radius, corresponding to 5 fm in the coordinate space, 
is used at all energies. The Z, A, excitation energy, angular momentum and velocity vector of 
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each cluster were calculated. Even after such a relatively long time, most clusters were in an 
excited state. In order to compare the simulated results to those of the experiments, the excited 
fragments were cooled down using the statistical decay code, GEMINI [10]. In this calculation, 
the C++ version of GEMINI was used. These events are referred to as the AMD + GEMINI 
events hereafter, whereas the events without the GEMINI calculation are called the primary 
events and referred to as the AMD events. The occurrence of the fusion reactions in the AMD 
+ GEMINI events is defined here by the emission of the fragments with Z > 6 in a given event. 

In FigJDthe calculated fusion cross sections, indicated by closed triangles, are compared to 
those of the experiments (open circles). The experimental data are taken from ref. [TTJ. The 
experimental data are reproduced well within the experimental errors above E cm > 3MeV in the 
absolute scale. The absolute cross sections predicted by the AMD simulations were calculated 
using the number of events generated in the given impact parameter range. At E cm < 3MeV the 
AMD simulation underestimated the fusion cross sections. In this energy range, the tunneling 
effect through the Coulomb barrier becomes important and in the present AMD formulation, 
this process is not incorporated. In the figure the formation cross sections of 19 F in the primary 
AMD events are also plotted by open square symbols. As discussed below, there are additional 
20 — 30% incomplete fusion contribution in the primary fusion process. 
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Figure 4. Fusion cross section for the 7 Li + 12 C. Circles represent experimental results and 
taken from [TT]. Squares are the primary of AMD results filtered by Z = 9. Secondary values 
of the AMD + GEMINI events filtered by Z > 6 are showed as triangles. 

In FigJSl the fusion channel distribution at the primary stage is shown as a probability 
distribution. Only the top three major channels are plotted. The 19 F formation and 15 N + a 
channel dominate the fusion reaction at all energies. Complete fusion occurs in about 80% of 
the cases at the lower incident energies and decreases to about 70% at the higher energies. The 
third channel contribution is from different reactions at different incident energies, but their 
probabilities are only a few % at most. 

In Figj6j the final exit channel distribution of the fusion reaction is plotted from the AMD 
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Figure 5. Primary major exit channel distribution at different incident energies. 
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Figure 6. Final exit channel distribution of the fusion reactions for the Li + 12 C. The blue 
histogram indicated the experimental values. The results of AMD+GEMINI are shown by green 
histograms, pn and pnn channels also include d and t decays. 



+ GEMINI events and compared with the experimental results of ref. [TT]. The relative cross 
section of the major decay channels is fairly well reproduced except the 15 iV + a channel. The 
suppression of this channel in the experimental results is not yet fully understood. Further study 
is now underway. 



3. SUMMARY 

The fusion cross section of the 7 Li + 12 C reaction was studied using the AMD and GEMINI 
codes. The AMD+GEMINI simulation reproduced the experimental total fusion cross sections 
reasonably well at E crn > 3MeV but underestimated it below that energy. The relative 
experimental exit channel distribution, except the 15 A + a channel, was well reproduced by 
the AMD+GEMINI simulation. 
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